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q . I. INTRODUCTION 

O 

There exist many papers dealing with algebraically special, expanding and twisting vacuum, pure radiation or 
Einstein-Maxwell fields. An extensive bibliography up to 1980 exists in [1]. Further results on Einstein-Maxwell fields 
may be found in [2] . Pure radiation fields have been studied extensively in [3-8] . Petrov type II fields are characterized 
by the non- vanishing of the second Weyl scalar ty 2 : 



* 2 = (m + iM) p 3 (1) 
Here to is the mass parameter, M is the NUT parameter and p is given by 



> , 

5: : p=-—^ ( 2 ) 

, where r is the coordinate along the null congruence of geodesies and £ is the twist [1], 

It has been noticed in different contexts that the condition M = simplifies the equations [2,4,5,7,9]. In the present 
t-H , paper we explore this condition with the help of the method proposed in [3], putting the emphasis on intrinsically 
time-dependent solutions. We clarify the general structure of the solutions obtained by separation of the traditional 
\ variables r, u, £, C where u is the retarded time and C, C span a two-dimensional surface. Known solutions are 
O^' systematized and new solutions are given. 

Twisting gravitational fields with M — generalize the classes of Robinson- Trautman (RT) [10] and Kerr-Schild 
bi). (KS) fields [11] and are physically realistic, their simplest representatives being the Schwarzschild, Kerr and Wadia 
' solutions. 

• i-h , In Sec. II the method of Stephani is applied to simplify the Einstein-Maxwell equations. In Sec. Ill the main 
field equation is reduced to a linear second order equation. Its general solution with separated variables is found and 
\ studied in Sec. IV. Sec. V is devoted to solutions linear in u, while in Sec. VI solutions with exponential behaviour 
- - ■ are discussed. Sec. VII contains some conclusions. 

II. THE METRIC AND FIELD EQUATIONS 

The standard form of the metric for expanding and twisting fields is [1]: 

ds 2 = - \dr + Wd( + Wd( + HQ] (3) 

ppP z L J 

fl = du + LdC, + Ld( 

It is determined by the r-independent real functions P, to, M and the complex function L: 

2iS = P 2 [pL - dl) (4) 

W = p^Lu + idT, (5) 
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H = -r (In P) u - (mr + MS - + y (6) 

# = 2P 2 i?e [9 (5 In P - L„)] (7) 
where <9 = — L<9„ and k is the Newton constant. The basic functions satisfy the following system of equations: 

(d- 3L U ) (m + iM) = -2 K p- 1 $?l>° (8) 

P- 3 M = ImddddV (9) 

^ = ~P [P-* (™ + iM)] u + P (ddddV) u - (ddV) u (ddV) u (10) 

(d-2L u )^ = (11) 

{d-L u ){p-^) + (p-^) u = (12) 
where = P, $J and $° are the essential parts of the Maxwell scalars 

$i = p 2 $? (13) 

$2 = + P 2 P (2L U - 5) $; + 2^p 3 P (SZ M - AS) $° (14) 
and n is either the energy density of pure radiation (then $i = $ 2 = ) or is given in the Einstein-Maxwell case by 

n 2 = 2k$°#° (15) 

Vacuum solutions have $i = $2 = n = 0. 

The form of the metric is preserved by certain coordinate transformations, one of which is the change of the retarded 
time 

u' = F(«,C,C) (16) 

P 1 = F-'P (17) 

L' = F U L - F c (18) 

E' = F^E (19) 

(m + iM)' = F~ 3 (to + iM) (20) 

Let us apply now to the system (8-12) the method of Stephani [3]. It works when either <j>° or vanishes. Then 
equation (8) is solved by 

m + iM = 4> 3 U (21) 
The complex field (f> is invariant under equation (16). It satisfies d<j> = 0, which stated differently gives 

L = T, 

When M = we can use equation (16) to transform to to a positive or negative constant to so that 
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<j> = m Q [u + iq ({,()] (23) 

L = iq c (24) 

Obviously L u = 0. This gauge differs from the most popular Kerr's gauge when P u = 0, but is very suitable when 
the NUT parameter vanishes. Eqs. (9) and (10) become 

ddddV = BBddV (25) 

n 2 

— = 3m P- 3 P u + PddddP - ddP ddP (26) 

with d = 8q — iq(d u . When P u ^ 0, n 2 can be made positive by the choice of to at least for some region of spacetime 
[1,4,9]. 
If $? = Eq. (12) gives 

$° = PG(0,C) (27) 
where G is an arbitrary function. If $2 = Eqs. (11) and (12) give 



$° = C(C)exp 



2i J (\nP) u q c d( 



(28) 



with C (£) being an arbitrary analytic function. This is a particular case of Theorem 26.2 from [1]. The expressions 
for the other functions in the metric simplify: 

S = P 2 Q (29) 

W = id (P 2 Q) (30) 

K = P 2 (Bd + dd) InP (31) 

H = -r (InP) u -(mor pp+ y (32) 

We have introduced the real, invariant under (16), function 

Q(C,C)=q C c (33) 

The Weyl scalars [12] contain a lot of terms with L u and also simplify in the gauge L u = 0. We quote only the leading 
terms: 

* 2 = m oP 3 (34) 
*3 = -p 2 P 3 dI + O (p 3 ) (35) 
* 4 = pP 2 Iu + O (p 2 ) (36) 

/ = P^ddP (37) 

What remains to be solved is the couple of equations (25) and (26) for P and q. For pure radiation Eq. (26) is just 
an inequality. 

When m + iM = (Petrov types III and N) Eq.(8) holds identically, but still a potential <j> may be introduced with 
the property dtp = and we can study the subclass of solutions satisfying Eq. (23) (with mo = 1). Then in all other 
equations we can set to = 0. 
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III. THE FIELD EQUATION (25) 



Eq. (25) is a linear equation of fourth order with respect to V. However, in the gauge L u = it becomes a linear 
equation of second order for P. This can be established with the help of the commutator 

[d,d\=2iQ8 u (38) 

Then Eq. (25) becomes 

Q cc -P + 2Q C -P C + 2Q C P- C + AQP a + 2i (q^Q c - q c Q^) P u + AiQ {q- Q P Cu - q c P- Cu ) + AQq c q- c P uu = (39) 

This equation characterizes the twisting solutions because when £ = it is trivial. It is non-linear in q except for 
time-independent solutions when only the first four terms remain. In this case Eq. (26) turns into the equation for 
type III RT solutions, containing only P. It is logical to solve first Eq. (26) for P and plug the result in Eq. (39) to 
find q. In the present paper time-dependent solutions will be discussed mainly. In this case Eq. (39) should be solved 
for P when q is given and the result placed in Eq. (26) to find n 2 for pure radiation solutions. 

Similar equations have been derived for time-independent fields in [1] (see Eq. (25.46)) and [2] (Eq. (5.12)). For 
time-dependent fields and concrete expressions for q such equations can be found in [7,9] . A linear equation of second 
order for (— 2E) X ^ 2 is found in [5] in terms of Cauchy-Riemann structures admitting Lie groups of symmetries and 
used in [4-6] to obtain pure radiation solutions. 

If P and q depend on £ and £ via a single function a (£, Q the terms in brackets in Eq. (39) cancel and it yields 

(Q aa a c a^ + Q a a C {) P + 4 (Q a a c a^ + Qa^) P a + AQa c a^ (P aa + q 2 a P uu ) = (40) 
Let ( — -j=(x + iy). Suppose that a = x and d y is a Killing vector. Then equation (40) becomes 

+ - A Q XX P = (41) 
and Q — \q xx . In the case of axial symmetry a = a = and Eq. (40) reads 

qlPuu + Pa, + ( ^P, + ^TtP = (42) 
cry 4erC l / 

Q = (43) 

The introduction of z = lncr as a variable greatly simplifies Eqs. (42) and (43) making them an analog of Eq. (41): 

qlQPuu + QPzz + QzPz + \QzzP = (44) 

Q = e- z q zz (45) 
Thus every solution with y-symmetry has an axisymmetric mirror with the same P and Q but as functions of z. 

IV. SEPARATION OF VARIABLES 

Let us search for twisting solutions by separating the variables P — p((,() f (u). Eq.(39) transforms into 

fun + Af u + Cf = (46) 



qcQ i(Q 





*■)} 




q-cv) . 



(47) 

C = (4g c g f Qp) _1 (Q^p + 2Q m + 2Q c p f + 4Qp u ) (48) 
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Eq. (46) is a well-known linear equation. It requires that A and C must be constant and possesses three types of 
solutions according to the sign of A 2 = 4C — A 2 : 

f = e' 4 ? (Cie^ + C 2 e~^ (49) 
/ = e-^(c lS in^ + C 2 cos^) (50) 

f = e-^(C lU + C 2 ) (51) 

where 1 = —A 2 > 0, A 2 > and A = respectively. Notice the common exponential factor if A ^ 0. It is quite 
interesting that damping exponential behaviour is generic for RT solutions, but it stems from the analog of equation 
(26) in the non-twisting case [13,14]. When Eq. (40) holds, A vanishes. There are still three types of solutions, one 
of them with exponential behaviour. When there is y-symmetry Eqs. (46) and (41) become 

fun +Cf = (52) 

Qpxx + QxPx+ (J^Qxx ~ CqlQ^J p = (53) 

In the case of axial symmetry equation (52) still holds, but equation (53) is replaced either by the same equation with 
x changed to z = In a or by 



p aa + (In oQ) a p a + 



4aQ 



p = (54) 



Let us discuss next Eq. (26) for two of the types of u -behaviour allowed by Eq. (39). Suppose first that P = pe cu . 
Introduce the function B — pe lcg . Then Eq. (26) reads 

n 2 = 6to c + 2p 2 (^e^ 4 ^) e 2lcq+4cu (55) 

V = BB ca - B c B- a (56) 

When mo = the solution is of type III only when ip ^ because dl = / B 2 . When mo 7^ and ip — the solution 
is of KS type because dl = 1^ = and / can be nullified by a coordinate transformation (' = g (() . This means 
ddP = (sec Eq. (37)), which can be lifted to ddV = 0. The latter is exactly the Kerr-Schild condition [1,3]. 
Suppose next that P = pu. When there is y-symmetry Eq. (26) yields 



n 



2 _ 6m 1 2 / 2 \ 4 r, 2 / ,^>\2 2 



u 2 

When there is axial symmetry Eq. (26) changes into 



■^P 1 (pPxxxx - P xx ) u* - 2p 2 (q x p x + Qp) u 2 (57) 



n 



, _ 6^0 + 2p2 [ 2mw + 4(TKW + a 2 ( ppaaaa _ p 2^ U A _ 2(J 2 [(g ff p2) J 2 u 2 (58) 



Finally, let us perform the separation of variables in the Maxwell equations. Eq. (28) simplifies 

$? = C(C)e X p[2i g (ln/)J (59) 
When f — u and / = e cu this becomes respectively 

= C (C) (60) 



= C (C) e 2lcq (61) 
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Eq. (12) or (27) separates as follows 

P- 1 ®" = N (C) e a{u+lq) (62) 

N (Q) being another arbitrary analytic function and a is a constant different from c in general. The energy density of 
a null Maxwell field (15) becomes 



k | JV(C) |" e zau (63) 



2P 2 



Obviously polynomial dependence of P in Eq. (26) is not allowed by (63). Consequently, null Maxwell fields can't 
induce solutions with P = pu. It is seen also from Eqs. (55), (57) and (58) that vacuum solutions of type II arc not 
possible because the different terms cannot cancel each other. 

V. SOLUTIONS LINEAR IN U 

These are solutions with C = in Eq. (46) or Eqs. (52), (53) and (54). Then these equations are identical to the 
equations for a u -independent P. Therefore we can use any time-independent solution of Eqs. (25) and (26) found 
in the literature which has L u = and axial or y-symmetry. We can also use time-independent solutions of Eq. (53) 
or Eq. (54) which are not solutions of Eq. (26). 

For example, let us find polynomial solutions with y-symmetry which capitalize on the only known vacuum RT 
solution for type III fields [1,10]. If we take L = ia; 6-1 , b ^ 1 then Eq. (53) becomes the Euler equation 

x 2 p xx + (b-2)x Px + ^(b-2)(b-3)p = (64) 

and its solutions are 

p=\x\ a (65) 

p = In | x | (66) 

p=\ x \^ [C*isin(^ln | x | ) + C 2 cos (fj, In | x Q] (67) 



when b < 3, b — 3 and b > 3 respectively. Here a = ^ (3 — b ± >/3 — b) and /i = \\fb — 3. The first of these solutions 
has been obtained essentially in [4] by a rather intricate procedure. Plugging Eq. (65) and q x = x h ~ x into Eq. (57) 
yields the inequality to be satisfied: 



L2 "'" ' 2 ( b - x± ^) +2a{a-i)(3-2a)u A -(2-V3~b) 2 \x\ 2b u 2 >Q (68) 

Throughout the paper we consider u > 0. The first term is always positive (we accept that mo > unless stated 
otherwise), the third is always negative, while the second changes sign and is positive for a < and 1 < a < 3/2. The 
third term has negative poles in x for any 6^0 and 6^—1 which cannot be compensated by the second term. Hence, 
pure radiation solutions of type III are not allowed and the first term should be present to compensate the negative 
poles with positive ones. This is achieved when b < 2. The plus sign in the formula for a must be taken in both cases. 
For growing u, however, the first term diminishes to zero. Thus the energy density is positive for b < 2 and small 
enough u. When b = the third term is x-independent but the second is negative, excluding again a type III solution. 
Type II solution still exists. Finally, when b = — 1 the third term vanishes but the second is still non-positive because 
a = 1 or a = 3. When a = 1 only the first term remains and n is x-independent, positive and bounded. This is a KS 
solution because ddV = —xq\j2 and can be turned into zero by a u-independent transformation of V. 

The second solution, given by Eq. (66), has no free parameters and leads to incurable negative poles and logarithmic 
singularities in n and is therefore unphysical. As a last comment , we may start with p = x b ^ 1 in Eq. (53) and obtain 
an Euler equation for Q with the corresponding three types of solutions. 

The simplest solution of Eq. (53) is probably when P is u-dependent only, i.e. p = po = const. Then Q xx = and 
P = pou, L = -j^x (c\x + 2c 2 ) where a are arbitrary constants and 
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2 6m 4 2 6 fcn \ 

n = 2p {) (ax + c 2 ) u (69) 

The energy density is positive for any x when c\ — and u is small enough. This is not a KS field. 

Another simple solution is obtained when p = e ax , a being some constant. Then Eq. (53) has constant coefficients 
and one of the solutions is q x = cxp (— 2ax). The energy density is positive, n 2 — 6m n /u and I = a 2 /2. Consequently 
the solution is of type II and is equivalent to a KS solution. 

Let us investigate next the axisymmetric solutions. As has been mentioned in the previous sections, Eq. (53) has a 
counterpart for axisymmetric fields with the same Q and P but with x replaced by z. However, q and n are different, 
being given by Eqs. (45) and (58) instead of q xx = 2Q and Eq. (57). The solution with p = p n has a mirror solution 

L = i( (c 3 lner + c 4 ) (70) 

n 2 = 5^0 _ 2c 2 p 4 u 6 

The energy density is positive for small u. When c 3 = it is positive everywhere and the solution is a KS field. 
Another solution has q — a, L = i(, Q = 1. Then Eq. (54) yields 

Paa + -Pa = (72) 

a 

and p = C5 In a + c 6 . It was found by a different method in [9] . When c 5 = it coincides with the previous solution 
with C3 = 0. 

A well-known solution is given by p = 1 + a/2. Eq. (54) becomes linear with respect to Q: 

(2 + a) oQ aa + (2 + 5a) Q a + 4Q = (73) 

and its solution is the hypergeometric function F (2, 2, 1, — a/2). It degenerates into a rational function for these 
values of its parameters and we have 

(74) 

(2 + af 

L = (75) 

(l + a/2) 2 

where k is an arbitrary constant, characterizing the magnitude of the twist (rotation). This is the Kramer solution 
[1,7,8,9] which is of KS type and has n 2 = 6mn/u. It represents a radiating Kerr metric. This is seen best in the 
gauge P u = 0. 

It must be stressed that any solution of the type P — pu can be transformed from the L u = gauge to the P u = 
gauge by choosing F = u 2 /2 in Eq. (16). Then P 1 = p, V = (2u) 1 ^ 2 L , m' = (2u) 3 ^ 2 m which transfers the 
M-dependence to the mass parameter and the twist. 

A generalization of the Kramer solution was undertaken in [7] for the same L and the ansatz P = A (a, u) (1 + a/2). 
It was shown that A satisfies a second order linear equation. The results of this paper guarantee that P also satisfies 
a second order linear equation. 

The mirror of the solution with p — e ax described above is p — e az and Q = k (1 — 2a) e~ 2az , a ^ 1/2. Hence 

P = a a u (76) 

L = ik(a- 2a (77) 

It is very similar to the Kramer solution for a = 1. For a = 1/2 , P is still given by Eq. (76) while L = ikQa^ 1 Inc. 
The energy density is 

n 2 = ^-25 h2a u 2 (78) 
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VI. SOLUTIONS WITH EXPONENTIAL BEHAVIOUR 



These are solutions with P = pe cu and C = — c 2 in Eq. (46) or Eqs. (52), (53) and (54). We shall give two solutions 
with y-symmetry. First we take q x = 2 (1 — b) x^ 1 , where the constant 6^1 and put this expression in Eq. (53): 



x 2 p xx + {b - 2) xp x + 



\{b-2) (6-3)+4c 2 (l-6) 2 x 2b 



p = (79) 



When 6 = this is the Euler equation, similar to Eq. (64) and has three types of solutions. The first two of them 
read 



p=|z| l±7 (80) 



p=\ x |5 In | x | (81) 

.1/2 



when c 2 < 3/16 or c 2 = 3/16 respectively. Here 7 = i (3 — 16c 2 ) . When b =/= the solution is given by Bessel 
functions: 

p = x 3 -^ [CiJ u (ex b ) + C 2 Y V (ex h )] (82) 



where v = V3 — b/2b, e = 2 | b — 1 | c/b. In all cases 

L = iV2{l-b)x b - 1 (83) 

These are essentially the solutions found in [4] where x-symmetry was used instead together with a special transfor- 
mation between a ^-independent and a u-dependent solution. 

The second solution to be presented is an analog of Eq. (69) i.e. P — e cu depends only on u. Then Eq. (53) may 
be integrated to 

1/2 \ 1/2 
Q = ^(b 2 --cV x + 2dq x j (84) 

where b > and d are constants. When d = Eq. (84) is solved by elliptic functions: 

L = i^cn(\x) (85) 
A 

where A = (|6 2 c 2 ) 1 ^ 4 and the modulus of cn is 1/2. It is clear from Eq. (84) that Q < 6/2. Eq. (55) yields 

n 2 = 6cm - 4c 2 Q 2 e 4cu (86) 

The second term is definitely negative and if m > 0, c > dominates for large u, no matter how big mo is. However, 
if m < 0, c < the first term remains positive while the second is damped exponentially for large retarded times. 
Choosing | mo | big enough we can arrange for n 2 > 0. 

Let us give next an example of an axisymmetric solution. Let us substitute in Eq. (54) q a = Q = a = const to 
obtain the generalization of Eq. (72): 



Paa + -Pa + Cl 2 C 2 p = (87) 
(7 



It is solved by Bessel functions e.g. 



P = e cu J (aca) (88) 
L = ia( (89) 

This solution was found and discussed in [9]. 

Finally, we give an example of a solution without y-symmetry or axial symmetry. We shall obtain the most general 
KS pure radiation field with P = pe cu . The results of Section IV and Eq. (56) show that 
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2ic B 



(90) 



where B satisfies tp = 0, namely 



P = (BB) 1/2 e c 



B 2 (lnB) cc - = S(C) 



(91) 



(92) 



S is an arbitrary analytic function. If S ^ we introduce the new variable to — j'' S (£') d(' and Eq. (92) transforms 
into the Liouville equation 



B 2 (InB) wC - = l 
Returning to the original variables, its general solution reads 

B = S (tf/* /(0+g(C) 

v//(0 C 5(c) c - 

The functions / and g are arbitrary. When 5* = the solution of Eq. (92) is 

B = /(C) 5(C) 



(93) 



(94) 



(95) 



It has zero twist i.e. it is a RT solution. Then we may set L = q = which means / (C) = g (C) and P = ff e cu . A 
coordinate transformation sets / = 1. Therefore the solution (95) is equivalent to B = 1. 

In [3,15] the general axisymmetric pure radiation KS fields were found. Interestingly enough, their w-dependence 
covers the three types (49), (50) and (51) with ^4 = 0. The exponential solution was generalized in [3] to a non- 
symmetric one with 

B = G(0+(K(() (96) 
G and K being arbitrary. Eq. (96) may be obtained from the general formula (94) by specializing to g (C) = C : 



B = f(0 



/(C), 



-1 1/2 



+ C 



s(0 

f(0c 



1/2 



(97) 



The general KS field with non-radiative Maxwell field has been given in [11]. It was remarked in [1] that no solution 
with a null Maxwell field has been found. Formulae (94) and (95) allow the study of this question in the context of 
solutions with separated variables. The energy density (63) tolerates exponential dependence. Eq. (55) reads in this 
case: 



3m c = kN (C) N (C) BB 



(98) 



Obviously B can't be given by Eq. (94) but a non-twisting solution with B = 1 is possible. Then P = e cu and 
$2 = (3m c/ k) 1 / 2 : . The mass parameter is constant, as usual, m — m . This is nothing but a special case of the 
Einstein-Maxwell solution of Robinson and Trautman [10] in disguise (see also Eq. (24.41) in [1]). The conclusion is 
that there is no twisting KS solution with a null Maxwell field and separated variables. 



VII. CONCLUSIONS 

In this paper we have carried out the idea of Stephani based on the introduction of an invariant potential for 
algebraically special, twisting and expanding gravitational fields. It works best in the important subclass of fields 
with vanishing NUT parameter. The system of Einstein-Maxwell equations is reduced tojthe couple of equations (25) 
and (26) for P and L where L is represented by a real function of two variables <z(C>C) (the L u = gauge). We 
have studied mainly pure radiational solutions for which Eq. (26) is an inequality. The main equation (25), originally 
of fourth order, becomes a linear second order equation for P (39), which further simplifies for the most common 
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symmetries. In some cases it is linear for q x too. Of course, it is not possible to enumerate all solutions of Eq. (39). 
The method of separation of variables was used to systematize most of the known solutions obtained in the past by 
a variety of different approaches. We have also investigated the region where the energy density of pure radiation is 
positive. Samples of new solutions were given to illustrate the application of Eqs. (39), (53), (54) and (55). Three 
types of behaviour with respect to the retarded time have been found: exponential, trigonometric and linear. Among 
the new solutions is the most general KS pure radiation field with exponential time behaviour. It was shown that the 
radiation field cannot be a null Maxwell field unless the twist vanishes. Finally, it should be mentioned that the case 
m = 0,M/0 is much more difficult because Eq. (9) remains non-linear in P. 
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